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Abstract 

Einstein-Gauss-Bonnet gravity (EGB) provides a natural higher dimensional 
and higher order curvature generalization of Einstein gravity It contains a new, 
presumably microscopic, length scale that should affect short distance properties of 
the dynamics, such as Choptuik scaling. We present the results of a numerical anal- 
ysis in generalized flat slice co-ordinates of self-gravitating massless scalar spherical 
collapse in five and six dimensional EGB gravity near the threshold of black hole 
formation. Remarkably, the behaviour is universal (i.e. independent of initial data) 
but qualitatively different in five and six dimensions. In five dimensions there is a 
minimum horizon radius, suggestive of a first order transition between black hole 
and dispersive initial data. In six dimensions no radius gap is evident. Instead, 
below the GB scale there is a change in the critical exponent and echoing period. 



1 Introduction 



Recent interest in string theory has popularized the study of higher dimensional and 
higher curvature gravity. The Einstein action has many desirable properties: it is second 
order in derivatives of the metric, ghost free when linearized about a flat background and 
obeys a Birkhoff theorem that yields a one parameter family of spherically symmetric 
black hole solutions. In four dimensions the Ricci scalar is the only curvature invariant 
with these properties but in dimensions greater than four it is possible to add higher 
order curvature terms in the form of so-called Lovelock polynomials [H El [3] . These terms 
contribute to the equations of motion while retaining no more than second derivatives of 
the metric. Moreover, they have been proven to be ghost-free [4] and obey a generalized 
Birkhoff theorem E] . 

The Lovelock action, /, written in terms of the Lovelock polynomials, C(p\ is given by 
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[n/2] 

d n x^f^g^a {p) C {p)) (1.1) 

p=0 



where K n = v / 8ttG^, G n is Newton's gravitational constant, g is the determinant of the 
metric, n is the number of spacetime dimensions, [n/2] refers to the largest integer less 
than or equal to n/2, SplZp* := 5^. ..5^ and 7t^ up(J is the Riemann curvature tensor. 
Q(p) are coupling constants of dimension (length) 2 (^ _1 ). The first two terms, £( ) and 
correspond to the cosmological constant and Einstein-Hilbert term, respectively, 
while C{2) is the Gauss-Bonnet (GB) term. It has been argued [7J that the GB term 
appears in the low-energy limit for strings propagating in curved spacetime. Here we 
focus on the simplest non-trivial theory, namely Einstein-Gauss-Bonnet (EGB) gravity, 
containing only the Einstein term and p = 2 GB term Q 

It has been known for quite some time that the spherically symmetric collapse of a 
massless scalar field minimally coupled to GR exhibits critical behaviour [8j . Specifically, 
numerical studies of black hole formation indicate that for any parameter in the initial 
data, A say, there exists a corresponding critical value A* such that for A > A* a 
black hole forms while for A < A* the matter disperses to infinity. Black holes with A 
just slightly bigger than A* are known as near critical black holes. Near criticality all 
geometrical quantities describing the black hole, such as its mass, obey a scaling relation 
of the form: 

ln(M BH ) = 7 \n(A - A*) + f(A - A*) (1.3) 

where / is a periodic function of its argument. The critical exponent, 7 and the period T 
of /, are universal in the sense that they are independent of the form of the initial data 



lr The addition of a cosmological constant should not affect the short distance behaviour that is the 
subject of our paper. 

2 7 depends only on the scaling dimension of the quantity considered. For example in 4D M scales 
as length so the radius and mass of the black hole have the same critical exponent. Curvature, on the 
other hand scales as (length) 2 , so its critical exponent is double that of mass. 
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or the specifics of the parameter A that is varied. This universality and the vanishing 



of Mbh a t criticality (cf. Eq.(1.3)) suggest a second order phase transition between 
the black hole and dispersive end states of the collapse. This fascinating behaviour was 
ultimately explained using renormalization group arguments in the context of radiation 
fluids[9j. For spherical massless scalar field collapse, the critical exponent and echoing 
period were obtained by Gundlach [10] from the properties of a critical, discretely self- 
similar zero mass black hole solution that behaves like an intermediate attractor in the 
space of solutions. 

The critical exponent 7 and period T do depend on the number of dimensions and 
the type of matter. In addition, the form of / can depend on the space-time slicing 
and the particular quantity that is being measured. For example, in Schwarzschild and 
null co-ordinates the radius of the horizon on formation yields an / that is well fit to 
a small amplitude sine wave. By comparison, in flat slice or Painleve-Gullstrand (PG) 
co-ordinates the periodic function that describes the scaling of the apparent horizon on 
formation exhibits large amplitude cusps [TU [15] . The difference can be understood by 
noting that in the former case one is measuring a quantity very close to the radius of 
the final event horizon, whereas in PG co-ordinates the apparent horizon is detected 
much earlier. The large amplitude cusps are likely due to that fact that near criticality 
in PG co-ordinates the apparent horizon forms at small radius and hence in the strong 
field region where such large fluctuations are expectecj^J On the other hand the scaling 
relation for the maximum value of the Ricci scalar at the origin for subcritical evolution 
is invariant and exhibits slicing independent small oscillations. 

The presence of a dimensionful constant in the action in general changes the above 
scenario, as verified for Yang-Mills collapse [IT], massive scalar field collapse [12J and 
massive gauge field collapse [13]. In massive scalar field collapse, for initial data whose 
width is smaller than the Compton wavelength of the scalar field, the usual second order 
phase transition is found, whereas in the other limit the phase transition exhibits a mass 
gap and is first order. It is clearly of interest to study the effects on Choptuik scaling of the 
Gauss-Bonnet parameter and higher order Lovelock coupling constants. Golod and Piran 
|20] recently presented such an analysis for the spherical collapse of massless scalar matter 
coupled to EGB gravity in five dimensions using double null co-ordinates. They found, 
as expected, that the Gauss-Bonnet term dominates the dynamics at short distances 
and destroys the discrete self-similarity characteristic of Choptuik scaling. Their work 
concentrated on the regime where the GB terms strongly dominated the dynamics. 

The purpose of the present work is to investigate further the critical collapse of a 
spherically symmetric, massless scalar field minimally coupled to five and six dimensional 
EGB gravity. We work in flat slice, or generalized Painleve-Gullstrand (PG), co-ordinates 
since they have several advantages over double null co-ordinates in the present context: 
They are regular at apparent horizons so that the simulations can run up to (and even 
past) horizon formation. Hence one can calculate the time and position of horizon forma- 
tion without having to stop the code at some arbitrary distance before horizon formation 
as would be necessary in Schwarzschild and null co-ordinates. More importantly, the 
cusp-like nature of the horizon scaling function in PG co-ordinates has the advantage 



3 We are grateful to Patrick Brady for suggesting this explanation. 
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of making the potential appearance of the periodicity in an equation such as (1.3) more 
obvious. It should be pointed out that one disadvantage of using PG co-ordinates is 
the lack of automatic spacial mesh refinement which occurs near horizon formation in 
null co-ordinates. As we will explain in the next section, the nature of the dynamical 
equations suggest that qualitative differences can occur in different numbers of space- 
time dimensions. It is for this reason that we investigate both five and six spacetime 
dimensions. 

We confirm some of the results in five dimensions [20], extend the analysis to six 
dimensions and obtain some surprising new results in both five and six dimensions. For 
all initial data and choice of parameter A that we examined there exists a critical value 
A* that separates black hole formation from dispersion. As expected, when the horizon 
forms far from the singularity the GR term dominates and the standard Choptuik critical 
scaling relation is found. Things change as one gets close enough to criticality to enter 
the region in which the GB terms dominate the dynamics. Near criticality the scalar 
field at the origin oscillates with a constant period T that converges as (A — A*) —> to 
a value that depends on the GB parameter as previously shown [20] . We find a different 
relationship between T and the GB parameter than in [20], albeit for smaller values of 
the GB parameter. 

In addition, we explore in detail the scaling in the GB dominated region. We find 
qualitatively different behaviour in five and six dimensions. In five dimensions there is 
evidence for a radius gap: in the supercritical region the radius of the apparent horizon 
on formation asymptotes to a constant value as criticality is approached from above. 
The maximum value of the Ricci scalar at the origin also appears to approach a constant 
value as criticality is approached from below. 

In six dimensions, the behaviour is qualitatively different. In the GB region the radius 



of the apparent horizon formation obeys a relationship similar to (1.3) but with different 
exponent and period. The maximum Ricci scalar also exhibits this same scaling relation 
with another scaling exponent, and small, but irregular oscillations. 

The rest of this work is organized as follows. In section [2] we describe the equations 
of motion which we derived using Hamiltonian formalism [18] . In section [3] we discuss 
the numerical implementation of the solution and describe the general methods used to 
obtain results. In section [| we give our results and conclude in section [5| 

2 Equations of Motion 

As stated above, we start with the action for a massless scalar field ip minimally coupled 
to the EGB action: 

1 = TT / d n x^~g (n + a {2) [U 2 - + Tl^U^] + k\ (V^) 2 ) , (2.1) 

ZK n J 

We use the ADM metric parametrization: 

ds 2 = -N 2 (x,t)dt 2 + A 2 (x,t)(dx + N r (x,t)dt) 2 + R 2 (x,t)dfl 2 (2.2) 
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where R is the areal radius. It is also useful to define the Misner-Sharp mass function [To], 
suitably generalized to EGBfTT]: 



M 



2kG 



2\ 2 



R n-z (! _ (DR) 2 ) + a {2) R n - 5 (l - (DR) 2 ) 



(2.3) 



where the n dimensional gravitational constant, G is defined as 2kG = 2k 2 / {n — 2)A n _ 2 
[Hj, Aj_ 2 is the surface area of an n — 2 dimensional sphere, k = 8(n — 3)/(n — 2) 2 , D is 
the 2 dimensional covariant derivative and (5(2) := {{n — 3)!/(n — 5)!)aj( 2 ). In the following 
we work in units in which 2G = 1. We work in flat slice co-ordinates x = R and A = 1 in 
which the equations of motion for the scalar field and its conjugate momentum are [TS] : 



ip = N 



IL 



R 



n-2 



+ 



N r 



(2.4) 



and 



n„ 



n I R"- V + ( ^ ) n„ 



(2.5) 



where dots and primes represent differentiation with respect to t and R respectively. 
Preservation of the gauge condition R = x in time determines the shift algebraically in 
terms of the lapse, the areal radius and the mass function: 



M := 



2kG 



R 



n—3 



N r 



This can be solved algebraically and yields, for EGB gravity: 



(2.6) 



~N 



\ 



R 2 



2a 2kGM 
a \ V + rV i?«- 3 ~ 



(2.7) 



where we have defined a := 2(n — 4)(n — 3)a(2) for purposes which will become obvious 
later. The sign of the inner square root in the above has been chosen to give the correct 
GR limit as a — » 0. 

The consistency condition A = determines the lapse via: 



R 



n-3 




NA & (Nr 

N \ R 2 \N 




(2.8) 



while the Hamiltonian constraint then takes the form: 



M' = 



TT 2 

2 KR^ 2 



R n - 2 ^' 2 



N r 
A" 



(2.9) 



Using (2.7) to replace N r /N by At in (2.9) provides a differential equation that can be 



solved for Ai and hence N r /N in terms of the scalar field and its conjugate momenta on 
each spacial slice. 
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Given the solution for Ai, one can look for apparent horizons by solving (2.3) for 
(DR) 2 = 0. In our coordinates, this becomes simply: 



AH := 1 



N r 
~N 







(2.10) 



where for ease of reference we refer to AH as the horizon function. For EGB, one can 



also use (2.6) and the above to obtain: 



M(Rah) = 



AkG 



(aR 



AH 



(2.11) 



Note that in 5D the first term is constant so that there is an algebraic lower bound on 
the black hole mass as the radius of the horizon goes to zero. 



Our goal is to solve the time evolution equations (2.4) and (2.5) for the scalar field 



and its conjugate momentum, with TV and N r /N on each time slice determined using 



(2. 7), (2. 8) and (2.9) and then use (2.11) to look for the formation of an apparent horizon. 



The actual time evolution equation as implemented in the code was obtained by ex- 



panding the derivative in (2.5) and replacing the derivatives of M! and N f using Eqs.(2.7) 



and (2.9). This gives: 



n,„ = n 




(2.12) 



Note that in five space-time dimensions the last term proportional to 1/R 3 in the square 
brackets above vanishes. One might therefore expect behaviour for n > 5 that is qual- 
itatively different from n — 5. It is for this reason that it is important to study higher 
dimensions. In the present paper we restrict consideration to five and six dimensions. 



3 Numerics and Methods 

The system is evolved using C++ code as follows: 

1. Initialize the spatial lattice. We set the lattice spacing to 10 -5 (unless otherwise 
stated) for the first 100 points near the origin and then slowly increase it to 10 -2 
at the 1200 t/l and final lattice point. 

2. Set up initial conditions. We initialized 11^ to zero and i/j to be either a Gaussian 
ipc or hyperbolic tangent ipn as follows 

ipc AR 2 exp 

where A, B and Rq are parameters. 



R — i?o 
B 



ipH A tanh 



R — Rq 
B 



(3.1) 
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3. At R — set N — 1, A4 — and use a subroutine to calculate N r /N using equation 



Eq.(2.7). Integrate N and Ai forward in R using equations Eqs.(2.7), (2.8) and 
(2.9). This is done using an RK4 method. Spatial derivatives are calculated using 
a central difference routine except at the boundaries where forward and backward 
difference are used. 



4. Integrate tj) and 11^ forward in time using equations (2.4), (2.5) and (2.7) employing 
an RK4 method. Stability is maintained by insisting that the size of the time step, 
At(t), is determined by 



fdR\ 



-l 



Af(f) < min j ^— j AR(R)}, (3.2) 

where AR(R) is the lattice spacing and ^ is the maximum value of either the 
ingoing or outgoing local speed of light. 

5. Monitor the apparent horizon function, AH := (DR) 2 = 1 — (^) 2 . At any point 
where AH = 0, there is an apparent horizon. When AH forms a minimum it 
signals that an apparent horizon is soon to form, so the time steps are diminished 
by a factor of ten. 

6. Calculate the Ricci scalar and the mass density. 

7. Repeat steps 3-6 until the formation of an apparent horizon or until the the field 
has dispersed. 

For comparison purposes it was important that the code could simulate collapse 
without the GB term, ie in the GR case. It is not possible to take this limit when 
numerically calculating N r /N using Eq.(|2.7D so an if statement was added to the routine 



which calculates N r /N in order to return N r /N = ^2kGA4/ R n ~ 3 when 5 = 0. 

When a is not zero a problem arises in the calculation of N r /N when AakG A4/ 'R n ~ l 
is sufficiently less than one. When this term is added to unity in the inner square root 



in (2.7), digits are lost and thus double precision can not be clai med. For this reason 
a 10 th order Taylor expansion of the inner square root in Eq.(2.7) was used in the case 



that 2a2kGA4/ R n ~ l < 0.1 Quadruple precision allowed for the investigation of overflow 
and underflow, as well as subtraction and addition round off errors. 

The code was capable of parallel processing, and many simulations were run on eight 
or more processors using the WestGrid and SHARCNET computing clusters. When 
generating data for mass and Ricci scaling plots the speed up was linear with the number 
of processors used, whereas for the binary search used to find critical values the speed 
up was logarithmic. 

The critical value of a parameter in the initial conditions is defined as the value of 
that parameter for which a black hole just barely forms. We first performed a binary 
search to find the critical value of A in Eq. fl3.ip . ^(t, R = 0) and A4 were then checked 
at late times confirm that they blew up for A slightly bigger than A* and remained finite 
for A slightly smaller than A*. We were able to get consistent results to 12 significant 
figures. The A* values for different values of a can be seen in figure [T] Interestingly 
the points are very well fit to straight lines. The above procedure, of course, also gives 
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£>* and which could also be varied. Using our values for A* we calculated the wave 
function at the origin as a function of PG time and used these plots to find the period 
of oscillation near criticality, as a function of the a. 



5D A* versus a 



6D A* versus a 
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(a) 5D (b) 6D 

Figure 1: A" as a function of a 

We simulated matter bounce and dispersal for 280 simulations (the number 280 cho- 
sen to optimize graph resolution and computing time) with A < A* and recorded the 
maximum value of the Ricci scalar 1Z at R = for each simulation. Plotting 1Z max as a 
function of A* — A with a log-log scale gives the Ricci scaling plot. Similarly we simluated 
collapse for 280 simulations with A > A* and recorded the radius of the initial appar- 
ent horizon, Rah- This procedure was repeated in five and six dimensions checking for 
scaling with both the A and B parameters in both the gaussian and tanh initial data of 
equation |3.1 1 to check for universality. Using gaussian initial data radius and Ricci scaling 
plots were created for a = 10 -8 , 10 _7 ,5 x 10 -7 , 10 -6 as well as the GR case a — to 
investigate the effects of the GB terms on the critical exponent, period and the existence 
of mass gap. 



4 Results 

4.1 Scalar Field Oscillations 

In general relativity the discrete self-similarity of the critical solution results in oscil- 
lations of the scalar field at the origin with ever decreasing period. The presence of 
the dimensionful Gauss Bonnet parameter breaks the scale invariance and the discrete 
self-similarity [20]. The scalar field oscillations at the origin near criticality approach a 
constant period that depends on the value of the GB parameter. Since it was difficult 
to get close enough to criticality to guarantee that the period had converged, we plotted 
the values as a function of log(cL4), where dA = \A — A*\. As seen in Figs |2(c)|3(c) 
the convergence was exponential and we used a best fit to determine the value of the 
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period T and its corresponding error for each value of a. The results are shown for 5 
and 6 dimensions in Figs |2(d)|3(d)[ Our results are qualitatively similar to those in [20] . 
namely 



T(n) oc a 1 



P 

(n) 



(4.1) 



Our exponents in five and six dimensions are: 



%) = 0.34 ±0.05 
/5 (6) = 0.24 ±0.08 



(4.2) 
(4.3) 



These both differ from the value of approximately 1/2 obtained in 5D by Golod and 
Piran [20] , who argued that f3 is one divided by the scaling dimension of the GB coupling 
coefficient. Intriguingly our results suggest a relationship of 



Note that the 6D plots show oscillations at late times which are likely due to the build 
up of numerical error. 

We note also that the range of a that we considered was between 5 x 10 -8 and 10 -6 , 
which is outside the range 4 x 10 -6 to 4 x 10 -4 considered by [20] . which may explain 
the discrepancy. We were restricted to smaller values of the GB parameter because our 
PG co-ordinate code did not allow us to get close enough to criticality for large values of 
a in order to reliably measure the period of the scalar field. 

4.2 Critical Exponents 

In GR there exist universal scaling relations whose properties are determined in part 
by the critical solution. We now present two different sets of scaling plots in the GB 
case. The first is the value of the logarithm of the apparent horizon radius Rah on 
formation as a function of log(oL4) as the critical parameter is approached from above 
(i.e. supercritical). The second is the log of the maximum value of the Ricci scalar 
at the origin as a function of log(cL4). We find as expected that if we are far enough 
from criticality that the curvatures stay small and the apparent horizon radius is large 
compared to the GB scale, we reproduce approximately the GR results: the curves are 
universal, with slope approximately equal to the GR critical exponent. The Ricci plots 
in this region are approximately straight lines with a small oscillation superimposed, 
whereas the radius plots show the large amplitude cusps observed in [T5l 02] . 

As the critical parameter is approached and we enter into the GB region things change. 
In the case of the Ricci plots, the GB region occurs when IZa > 1, whereas for the radius 
plots it can be defined by the simpler relation Rah < V&- The boundary between the 
two regions is indicated in all the scaling plots by a horizontal dashed line. 

The radius plots Figs. ^ and ^ continue to exhibit cusps, but with a decreased 
period and slope. The Ricci plots are also similar in the GB region to the GR region 
in that they are approximately straight lines with oscillations superimposed. However, 
the slope changes quite suddenly when the transition from GR to GB is made. The first 



/3 (n) = l/(n - 2) 



(4.4) 
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Figure 2: Scalar field oscillations 



important point is that the scaling plots are universal even in the GB region. This is 
illustrated for both 5D and 6D in Fig.Q. There are qualitative differences in the scaling 
plots between 5D and 6D so we will now discuss the two cases separately. 

In the case of 5D there is evidence that the slope of the radius plot decreases con- 
tinuously until a minimum radius is reached, i.e. that there is a radius gap. This is 
most evident in Fig.( |6(d)j ) but also appears to be the case in Q6(e)D . In the remaining; 5D 
figures the numerics did not allow us to probe deeply enough into the GB region to fully 
observe this. 

A radius gap is not unexpected given the presence of the dimensionful GB parameter. 
Note that we focus on a radius gap instead of a mass gap because in 5D the former is 
trivial in light of (2.11). 

The Ricci plots, initially approximately straight, change slope quite suddenly as one 
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(b) 6D,^(0,t) near criticality, a = 10 6 . 
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(d) Period of ip(0,t) in 6D near criticality 
as a function of GB parameter 



Figure 3: Scalar field oscillations 



moves from the GR to the GB region, and then remain constant over a small range of 
log(oL4). The slopes are given in Table [l] As criticality is approached the slope of the 
Ricci plot gradually decreases, suggesting that there is a maximum value to the Ricci 
scalar at the origin. This differs from GR, in which the critical solution is singular and 
the Ricci scalar at the origin increases indefinitely as criticality is approached. 
We emphasize again that these features are universal. 

In 6D things are different. There is no evidence of a radius gap in the radius scaling 
plots, and the slope of the Ricci plots remains constant until we reach the limits of 
numerical accuracy. Thus it appears that there is a transition to a new set of scaling 
exponents, which are plotted in Table ([2]). Note that numerical uncertainties make the 
first and last entries in each column unreliable. The exponents are different for Ricci 
and radius scaling, but the absolute value of both appear to increase with decreasing a. 
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a 


5D Ricci Scaling 


10" 6 


-1.426 ±0.074 


5 x 10" 7 


-1.573 ±0.076 


10- 7 


-1.577 ±0.028 


10" 8 


-3.397 ±0.049 



Table 1: 5D Ricci scaling exponents in GB region. 



a 


6D Ricci Scaling 


6D Radius Scaling 


10" 4 


-1.488 ±0.128 


0.257 ± 0.002 


10" 5 


-1.433 ±0.016 


0.207 ±0.002 


10" 6 


-1.619 ±0.021 


0.313 ±0.002 


5 x 10" 7 


-1.814 ±0.016 


0.476 ± 0.002 


lO" 7 


-2.029 ±0.027 


0.417 ±0.002 



Table 2: 6D Ricci and AH radius scaling exponents in GB region. 



Moreover a log-log plot of the three reliable Radius vs Ricci exponents (Fig. [9]) reveals 
that they are related by: 

1{mcci) « -(2.24 ± 0.04) x 7 ^ S ) 2 (4-5) 

This is to be compared to the GR case in which the relation is determined purely by the 
dimension of the two quantities: 

l(Ricci) = -^(Radius) (4.6) 
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Radius Plot for a=5e-07 in 5 Dimensions 
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(b) 5D a = 5 x 10 T , amplitude and width 
shifted to lie on top of each other 
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(c) 5D a = 5 x 10 T , Radius Plots Super- 
imposed 
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(d) 5D a = 5 x 10 7 , Ricci Plots Superim- 
posed 



Radius Plot Showing Universal Scaling 6D 



Ricci Plot Showing Universal Scaling 6D 




i i i 1 1 i i i i 1 1 i 1 1 1 1 i i 1 1 i i i 1 1 i i 1 1 i i i 1 1 i i i i i 1 1 i i i 
-11 -10 -9 -8 -7 -6 -5 -4 -3 -2 
log io (|par o -par*|) 

(e) 6D, a = 10 -5 , Radius Plots Superim- 
posed 
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(f) 6D, a = 10 5 , Ricci Plots Superimposed 



Figure 4: Universality in 5 and 6D 



Radius Plot for a=0e+00 in 5 Dimensions 



Radius Plot for ct=1e-08 in 5 Dimensions 
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(a) GR, slope= 0.413 

Radius Plot for a=1e-07 in 5 Dimensions 
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Radius Plot for ot=5e-07 in 5 Dimensions 
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(d) a = 5 x 10- 7 
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(e) a = 1(T 6 

Figure 5: Radius Scaling Plots - 5D. The lines represent the best-fit tangents to the 
curves in their respective regimes. 
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Ricci Plot for a=0e+00 in 5 Dimensions 



Ricci Plot for a=1e-08 in 5 Dimensions 





(a) GR, slope= 0.826 

Ricci Plot for a=1e-07 in 5 Dimensions 
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Figure 6: Ricci Scaling Plots - 5D. The lines represent the best-fit mean slopes of the 
curves in their respective regimes. 



15 



Radius Plot for a=0 (GR) in 6 Dimensions 



Radius Plot for a=1e-07 in 6 Dimensions 
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(a) GR, slope= 0.43 

Radius Plot for a=5e-07 in 6 Dimensions 
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Figure 7: Radius Scaling Plots - 6D. The lines represent the best-fit tangents to the 
curves in their respective regimes. 
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Ricci Plot for a=0 (GR) in 6 Dimensions 



Ricci Plot for a=1e-07 in 6 Dimensions 




(a) GR, slope= 0.43 

Ricci Plot for a=5e-07 in 6 Dimensions 
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Figure 8: Ricci Scaling Plots - 6D. The lines represent the best-fit mean slopes of the 
curves in their respective regimes. 
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Critical Exponents GB Region 6D 




Figure 9: Plot of Radius exponents vs Ricci 6D 
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5 Conclusion 



We studied the effects of the GB term on the dynamics of the collapse of a massless scalar 
field minimally coupled to gravity in five and six spacetime dimensions. The GB term 
destroys the self-similar behaviour, as demonstrated by the fact that near criticality the 
scalar field at the origin oscillates with a constant period. The period in five dimensions 
is proportional to roughly the cube root of the GB parameter and as the fourth root 
in six dimensions. While the 5D results differ from those in [20] it must be emphasized 
that we have explored a different range of GB parameter, and this may account for the 
difference. 

We also showed the existence of modified, but still universal, horizon and Ricci scaling 
plots near criticality. We found evidence for the existence of a radius gap in five dimen- 
sions but not in six dimensions. This qualitative difference is not completely unexpected. 
As mentioned below Eq.(2.12), the time evolution equation in five dimensions is special, 
containing one less term than in the higher-dimensional cases. It may also be useful to 
note that qualitative differences exist between five and six dimensions with regard to the 
stability of black holes under gravitational perturbations j22j [23]. Small five dimensional 
GB black holes are unstable with respect to scalar gravitational perturbations, whereas 
in six dimensions it is the tensor mode that yields an instability. 

It is clearly of interest to confirm our results with further simulations and to try to 
understand analytically the source of the new scaling behaviour. 
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